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We propose an unexplored class of absorbing materials for high-efficiency solar cells: heterostruc- 
tures of transition-metal oxides. In particular, LaV03 grown on SrTi03 has a direct band gap 
~ 1.1 eV in the optimal range as well as an internal potential gradient, which can greatly help 
to separate the photogenerated electron-hole pairs. Furthermore, oxide heterostructures afford the 
flexibility to combine LaV03 with other materials such as LaFe03 in order to achieve even higher 
efficiencies with band-gap graded solar cells. We use density-functional theory to demonstrate these 
features. 

PACS numbers: 88.40.fh, 73.20.At, 78.20.Bh 



The stunning discovery by Ohtomo and Hwang [1, 2] 
that the interface between the band insulators SrTi03 
and LaAlOg (or LaTiOg) can become conductive above 
a critical thickness of the latter opened the research field 
of layered oxide heterostructures. The interest comes not 
only from novel physical effects that are absent in the 
constituent bulk materials, but also from the perspective 
of designing and tuning specific properties to achieve de- 
sired functionalities. This flexibility becomes even more 
pronounced when materials with partially filled d-shells, 
and hence strong electronic correlations, are involved 
[3, 4]. In particular, these novel oxide heterostructures 
raise hope to surpass state-of-the-art semiconductors for 
specific applications. In this Letter, we propose high- 
efficiency solar cells as such an application. 

The performance of photovoltaic systems is character- 
ized by the efficiency t] of conversion of the incident pho- 
ton energy to electrical power. A natural limit for 77 is set 
by the optical gap A of the absorbing material: photons 
with energy hio < A are not absorbed at all; those with 
huj > A may be absorbed but each only contributes an 
amount of energy equal to A, with the difference being 
lost to relaxation processes, i.e., heat generation. This 
consideration, combined with further loss channels, leads 
to the famous Shockley-Queisser limit [5], which gives an 
upper bound for ry as a function of A with rj < 34% in 
the optimum range A = 1-1. 5eV. In the case of semicon- 
ductors, a great effort now concentrates on the reduction 
of so-called Shockley-Read-Hall electron-hole recombina- 
tion caused by defect-induced "trap" states and, in gen- 
eral, any intrinsic recombination mechanism. Even small 
improvements on the final efficiency of a solar cell have 
an immense economical and environmental impact, and 
are intensively pursued by industry [6, 7]. 

Why oxide heterostructures ? There are four prop- 
erties that make layered oxide heterostructures, and in 
particular LaVOglSrTiOg, a promising candidate for ef- 
ficient solar cells (cf. Fig. 1): (i) an intrinsic electric field 




Figure 1: (Color online) Schematic of a gap-graded oxide- 
heterostructure solar cell. The cell shown here combines two 
absorber materials, generically LaM'Os and LaM03 (M and 
A/' being different transition metals); the compound with 
the larger band gap A' > A should face the sun. Incom- 
ing photons with energy huj > A' can be absorbed in the top 
material and contribute an energy A', whereas those with 
A < fioj < A' traverse the first part, which is transparent to 
them, and can be absorbed in the lower material, contribut- 
ing an energy A. Once an electron-hole pair is dissociated, 
the built-in effective electric field E* will drive the carriers to 
their respective contacts, which are naturally provided in the 
form of conducting surface and interface layers. 



emerges in the photoabsorbing region, which may effi- 
ciently separate photoexcited electrons and holes; (ii) the 
band gap of 1.1 eV is direct and in the ideal energy range 
for harvesting the sunlight reaching the Earth's surface; 
{Hi) interfaces and surfaces can become metallic on a 
thickness of about one unit cell, and those layers natu- 
rally allow for extracting the charge carriers; and (iv) one 
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Figure 2: (Color online) multi-layer LaV03|SrTi03 het- 
erostructure with 4 VO2 planes. Superimposed on the 
relaxed structure, the red/blue (light/dark gray) lobes 
show isosurfaces of the positive/negative part of the 
spin-density n-f — nj,. Owing to the periodic repetition 
of the cell, both an n-type (left) and a p-type (right) 
interface appear. The substrate is only partly shown; 
altogether there are six Ti layers. As in bulk LaV03, 
the V-d electrons show AF-C spin order; for the central 
V layers, the AF-G orbital order is similarly preserved. 



can flexibly combine different materials [8-12] and grow 
heterostructures with multiple band gaps, e.g., alternate 
LaV03 with LaFe03 in order to realize a so-called band- 
gap graded design [13]. By means of density-functional 
theory (dft) calculations we establish these properties 
of LaVOglSrTiOg and, given the enormous flexibility of- 
fered by oxide heterostructures, we suggest them as a 
new candidate for efficient photovoltaic devices. 

LaAlOglSrTiOg can be considered the prototype for 
the kind of heterostructure we consider here, and has 
attracted by far the most attention [14-18]. However, 
for our purposes LaVOgjSrTiOg is more suitable since 
LaAlOgjSrTiOg has a band gap of several electron volts, 
much too large for photovoltaics. Both heterostructures 
share a feature which is most interesting in relation to 
the electron-hole recombination problem: a polar inter- 
face with an intrinsic electric field. In LaVOg both La 
and V have nominal valence -|-3 such that the LaO planes 
are positively charged and the VOj planes are negatively 
charged, inducing an electric field between these layers 
(depicted schematically in Fig. 1; for quantitative results 
see Fig. 4 below). In the substrate SrTiOg, on the other 
hand, Sr has nominal valence +2 and Ti has +4, such 
that the SrO and TiOj planes are both charge neutral. 
Based on these simple considerations one expects a po- 
tential gradient induced by the polar discontinuity be- 
tween LaVOg (or LaAlOg) and SrTiOg [19]. dft calcula- 
tions for LaAlOglSrTiOg show that, in order to avoid this 
so-called polarization catastrophe, the polar discontinuity 
is partially compensated by electrons transferred from 
the top of the heterostructure to the interface, whereas 
a substantial gradient persists. Recently, a finite poten- 
tial gradient has also been measured experimentally [20] . 
That said, such measurements are challenging, in partic- 
ular because of the difficulty of controlling oxygen vacan- 
cies in LaAlOg (or, equally, LaVOg), which may greatly 
influence the gradient. 

The built-in electric field will support the separation 
of photogenerated electron-hole pairs and will drive the 
dissociated charge carriers selectively to the respective 
contacts for extraction. LaVOglSrTiOg combines this 
feature with a near-optimum direct band gap of 1.1 eV. 
Therefore, we concentrate on this structure in the present 
Letter. 
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Figure 3: (Color online) Density of states of a multi-layer 
LaV03|SrTi03 heterostructure. Each of the four VOj lay- 
ers (V1-V4) contains two inequivalent V atoms (Va and Vb) 
carrying opposite spin; the majority-spin V contribution is 
shown for each atom. For the Ti contribution (we show only 
the interfacial layer), the two spins as well as the two sites 
are essentially identical. These are the only relevant contri- 
butions around the Fermi level. Note the layer-by-layer shift 
of the bands which indicates the potential gradient, and the 
appearance of conducting states (of d^y character, only in one 
spin channel) on the interfacial V. 

Another advantage of the proposed oxide heterostruc- 
tures for solar-cell applications is that electrical contacts 
to collect the photogenerated charge carriers are natu- 
rally provided. The ("n-type"; see below) LaVOglSrTiOg 
interface itself is metallic; in fact, it would be difficult to 
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Figure 4: (Color online) Potential gra- 
dient in multi- layer LaV03|SrTi03 het- 
erostructures. The curves (■/ jk distin- 
guish structure size) track the energy of an 
O-ls state through the TiOj (layer < 0) 
and VO2 (layer > 0) layers, providing a 
measure of the potential gradient. Above 
the data, the layers, interfaces (IF), and 
effective electron potential V* are shown 
schematically. Periodicity forces the poten- 
tial to return to zero in the SrTi03 part 
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contact it otherwise. Our dft calculations show that the 
surface layer is also metallic and, hence, suitable to ex- 
tract the holes. However, surface defects and edges might 
localize the holes. Hence, careful surface preparation or 
even an additional metallic layer such as SrV03 or other 
metallic contacts will be necessary on the surface side. 
Ah initio simulations of oxide solar cells We 
study structures with different numbers of LaO and VO2 
layers grown on the (001) surface of SrTi03 consid- 
ering both (i) periodically repeated arrangements con- 
taining one "n-type" (TiOjl [LaO]''') and one "p-type" 
(SrO|[V02] ) interface and (ii) an n-type (Ti02|LaO) 
interface with a slab of vacuum on top of the LaVOg 
part to break periodicity. These setups are termed (i) 
"multilayer" and {ii) "thin-film" geometries. 

The DFT-I-U calculations were performed using the full- 
potential linearized augmented plane-wave code WiEN2k 
[21], with the PBEsol [22] and mBJ [23] exchange- 
correlation potentials and a local Coulomb interaction 
term U on Ti, V and Fe. The absorption coefficients in 
Fig. 5 were calculated using the OPTIC code for WiEN2k 
[24]. See Supplemental Material [25] for technical details. 

We begin by showing one representative 
LaVOgJSrTiOg structure, including the results from 
structural optimization and electronic self-consistency, 
in Fig. 2. As in the bulk [26, 27], the Y-d electrons show 
staggered spin and orbital order. The spin order is of 
the C-type (ferromagnetically coupled antiferromagnetic 
planes), while the orbital order is of the G-type (alter- 
nating in all directions) [28] . The V and Ti contributions 
to the density of states for this case are shown in Fig. 3, 
in the energy range of the V-d bands. These bands shift 
layer by layer, a first sign of the potential gradient. The 
gradient can be evidenced more clearly by tracking the 
energy of a core level throughout all layers, as shown in 
Fig. 4. 

We estimate the potential slope in the LaVOg region 
to be 0.08 eV/A (see Fig. 4). A similar potential slope 
is observed in the LaVOg region of thin-film structures 
(see Supplemental Fig. 1 [25]). As the photocarrier exci- 
tation takes place in the LaVOg region, such a potential 
gradient inside LaVOg regions helps to efficiently sepa- 



rate photoexcited electrons and holes. 

For the multi-layer heterostructure with four V lay- 
ers, we observe the appearance of states at the Fermi 
level at the interfaces between LaVOg and SrTiOg. This 
confirms the scenario of metallic interfaces due to the 
electronic reconstruction found also in LaAlOg] SrTiOg 
within DPT. Yet, contrary to LaAlOg] SrTiOg, the carri- 
ers reside mostly on the V rather than the Ti, in accor- 
dance with recent experimental results [29]. The criti- 
cal thickness that we find for multi-layer LaVOg [SrTiOg 
heterostructures at the DFT-I-U level is four. This result 
compares well with what has been reported in one ex- 
periment [30] and it is smaller than what was found in 
another experiment [29]. While any solar cell must have 
a band gap in order to generate electrical energy, con- 
ducting states that stay confined to the interface may in 
fact prove useful in extracting the photogenerated charge 
carriers. 

Let us now turn to a central quantity for solar cell 
applications, the optical absorption (i.e. a in the Beer- 
Lambert law I{r) ~ e""*", see Supplemental Material 
for details [25]), shown in Fig. 5. A major advantage 
of LaVOg I SrTiOg heterostructures over the current stan- 
dard absorber material is that, contrary to Si, the band 
gap in LaVOg is direct, such that photons carrying the 
band-gap energy can create electron-hole pairs without 
the aid of phonons or another indirect scattering pro- 
cesses. Furthermore, across most of the solar spectrum, 
bulk LaVOg compares favorably with CdTe, a direct-gap 
material currently used for high-efficiency, thin-film solar 
cells. 

The validity of the absorption coefficient computed 
within DPT is confirmed by a comparison to experimental 
data on bulk LaVOg [27]. Note that the first sharp peak 
in the theoretical curve is very sensitive to details of the 
crystal structure (in the limit of a cubic unit cell, it van- 
ishes by symmetry). This might explain why it is absent 
in the experiment and in the heterostructure calculation. 

LaVOg I SrTiOg heterostructures have similar optical 
properties as bulk LaVOg, with the additional freedom 
that individual layers can have different band gaps owing 
to layer-dependent distortion and resulting crystal-field 
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Figure 5: (Color online) Absorption coefficients and solar 
spectrum. The lines show the absorption coefficients of bulk 
LaV03 and LaV03|SrTi03 compared to experimental data 
on bulk LaVOg [27] and a calculation for CdTe, which is 
widely used in current high-efficiency solar cells. In the back- 
ground, the solar spectrum as measured on the Earth's surface 
is shown (standard global air mass 1.5, in arbitrary units). 



splittings. Another advantage of oxide heterostructures 
is that we can flexibly combine LaV03 with a second 
transition metal oxide with a somewhat larger band gap, 
such as LaFeOg, improving the conversion efficiency in 
the high-energy region of the solar spectrum. Thus one 
can construct gap-graded structures, in which the lay- 
ers facing the sun have larger gaps than the layers un- 
derneath. The SrTi03 substrate can be on either side 
because it is transparent to visible light. This scheme 
efficiently reduces the losses due to the magnitude of 
the band gap and provides the primary means to exceed 
the Shockley-Queisser limit. We perform dft-|-u calcu- 



lations for such a LaFeOg jLaVOgjSrTiOg heterostructure 
and find a layer-dependent gap of 2.2 eV for the LaFeOg 
and 1.1 eV for the LaVOg part, in close agreement with 
the bulk values. See Supplemental Fig. 3 for details [25]. 
Note added in proof Manousakis [Phys. Rev. B 
82, 125109 (2010)] reports that, in a Mott insulator, a 
single photon may create multiple electron-hole pairs by 
impact ionization. This effect could further increase the 
efficiency of the LaVOgjSrTiOg solar cell proposed here. 
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The DFT calculations were performed using the 
full-potential linearized augmented plane-wave code 
WlEN2k [21], with the PBEsol exchange-correlation po- 
tential [22] and a local Coulomb interaction term U on 
Ti, V and Fc {Uti = 9.8 cV, Uy = 3cV, Ufc = 5cV, i.e. a 
GGA-I-L'' scheme. The U values have been selected to re- 
produce the experimental band gaps in the bulk, which 
is essential in this case since the heterostructures would 
become conducting too soon (i.e. the critical thickness 
would be too small) if the gaps were smaller. 

This is also the rationale for the large C/ti: It is neces- 
sary to shift the Ti-d states up far enough; but because 
they are almost empty, this drastic value is needed to 
achieve the desired effect. This is analogous to the case 
of LaTiOg, where a large U on the empty La-/ states is 
necessary to shift them away from Ti-d. 

An independent verification of the results achieved 
with this 17x1 comes in the form of a calculation using the 
parameter-free modified Becke-Johnson (mBj) exchange- 
correlation potential [23], which we find to qualitatively 
agree with the corresponding DFT-|-u calculation. 

Formally, periodic boundary conditions are in force in 
all our calculations; however, in the thin-film case the 
vacuum acts to separate LaVOg and SrTi03, effectively 
imposing open boundary conditions. 

The lattice constants of SrTiOg, LaVOg, and LaFeOg 
arc quite similar when one takes into account that the 
unit cells of the latter two are enlarged with respect to 
the primitive cubic perovskite cell due to the reduced 
symmetry; they are doubled in the z-direction, and form 
a v^-cell in the xy-plane. We take this setting for 
the heterostructure unit cell in the xy-pl&ne, thus there 
are two Ti/V/Fe atoms per plane. We fix the in-plane 
lattice parameter to the value in the SrTiOg substrate 
(this corresponds to the experimental situation), while 
the extent in z-direction for the multi-layer structures 
corresponds to the sum of the bulk layer separations. 
For the corresponding thin-film structure, this value is 
enlarged by 20 Bohr 10.6 A. Within this volume, all 
atomic positions are relaxed. For the numerical values, 
see Supplemental Table I. 

In the bulk, both LaVOg and LaFeO;^ show antifer- 
romagnetic (AF) spin and orbital order at low temper- 
atures: LaVOg has AF-C spin and AF-G orbital order 
while LaFeOg has AF-G order in both channels; dft 
predicts this in accordance with experiment. This raises 



Supplemental Table I: Lattice constants of bulk SrTiOg (cu- 
bic), LaVOj (orthorhombic), and LaFeOg (orthorhombic) . a 
and b arc the in-plane y/2 lattice parameters, while c is the 
BOj-plane separation. 





a/ A 


fe/A 


c/A 


SrTiOg 


5.523 


5.523 


3.905 


LaVOg 


5.539 


5.560 


3.907 


LaFeOg 


5.554 


5.568 


3.928 



the question which spin order prevails in the heterostruc- 
tures. For LaVOg] SrTiOg, we focus on AF-C both in 
analogy to the bulk, and because it has the lowest en- 
ergy in DFT (among AF-A, AF-C, AF-G and FM) in two 
test cases. 

When a second magnetic material comes into play with 
LaFeOg, the situation becomes more complicated and 
three different spin orders have to be considered in gen- 
eral: within LaVOg, within LaFeOg and at the interface. 
In the LaFeOglLaVOgjSrTiOg calculation presented here, 
there is only one layer of each material, therefore we 
need to consider only the coupling between LaFeOg and 
LaVOg, which we choose antiferromagnetic in the inter- 
est of comparability with the LaVOg-only calculations, 
resulting in "AF-C" spin order across V and Fe. 

The absorption coefficient for LaVO-j was calculated 
within GGA+U as described above, using the program of 
Ref. 24 while for the semiconductors the mBj potential 
was used [23] , which yields band gaps in good agreement 
with the experimental values in these cases. The calcu- 
lated absorption coefficients have been checked regarding 
convergence in the density of fc-points used to sample the 
Brillouin zone. For the heterostructure, up to 10 x 10 x 1 
/c-points were used. 

The experimental absorption coefficient for LaVOg was 
calculated from the reflectivity measured in Ref. 27 using 
Kramers-Kronig relations following Ref. 31. 
Gradient in the thin-film case. In analogy to Fig. 4 
in the main paper, we show the behavior of the 0-ls level 
in thin- film heterostructures in Supplemental Fig. 1. 
Absorption coefficients. As stated in the main text, 
the absorption coefficient a(w) is defined via the expo- 
nential dampening of the light intensity as a function of 
distance travelled through the material. However, in ad- 
dition to the frequency dependency, the absorption also 
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Supplemental Figure 1: Potential gradient 
in thin- film LaV03|SrTi03 heterostruc- 
tures. The curves (■/•/! distinguish 
structure size) track the energy of an O- 
Is state through the TiOj (layer < 0) and 
VOj (layer > 0) layers, which provides a 
measure of the potential gradient. Above 
the data, the layers, interfaces (IF), and 
effective electron potential V* are shown 
schematically. The potential in the SrTiOg 
part is almost flat, showing only a slight 
buckling which we attribute to band bend- 
ing at the interface. 
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Supplemental Figure 2: Contributions to absorption coeffi- 
cients. Both in bulk LaVOg and the heterostructure, a signif- 
icant asymmetry is seen between the xx and the zz compo- 
nents, while yy is similar to xx. The off-diagonal components 
are rather small in comparison with the diagonal ones. 



depends on the light's polarization, leading to a matrix 
aij analogous to the dielectric tensor. In cubic crystals 
(such as CdTe and GaAs), symmetry reduces the tensor 
to only one independent component, an, and this is what 
is plotted for these materials in Fig. 5 in the main paper. 

On the contrary, in LaVOg and LaV03|SrTi03, the re- 
duced symmetry allows six independent entries. Which 
of these may contribute to the actual absorption de- 
pends of course on the light's angle of incidence with 
respect to the crystal. Since the heterostructure is grown 
along the z-direction, in the first approximation the po- 
larization will be in the xy plane. However, because 
some of the light will come from other directions due to 
scattering, and because heterostructures are often grown 
on a stepped substrate, some z-component will also be 
present, and all may contribute. Therefore, we plot- 
ted a suitably averaged quantity a := |(X]i '^ij) 
in Fig. 5. Supplemental Fig. 2 shows some of the indi- 
vidual contributions. The optical absorption for the het- 
erostructure was calculated in the multi-layer geometry 
with two VO2 layers. 

Band-gap grading with LaFeOg. In this section. 
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Supplemental Figure 3: The relevant contributions to the 
density of states around the Fermi level for the "gap-graded" 
LaFe03|LaV03|SrTi03 structure. In this case, the two V/Fe 
sites give almost identical contributions (in opposite spin 
channels). A p- and an n-type interface appear, as marked, 
due to periodic boundary conditions. 



we give details on the calculation combining LaV03 and 
LaFe03 referenced in the main text. The structure con- 
tains 3 Ti, 1 V and 1 Fe layer, with two metal atoms per 
layer. The projected density of states for Fe, V and the 
Ti layer at the n-type interface is shown in Supplemental 
Fig. 3. The differing magnitudes of the Fe-0 and V-V 
gaps are clearly seen. 

It should be emphasized that LaV03 -I- LaFe03 is only 
one example for many different combinations that may 
be envisaged within the flexible framework provided by 
oxide heterostructures. 
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